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We consider non-interacting bosonic excitations in disordered systems, emphasising generic features 
of quadratic Hamiltonians in the absence of Goldstone modes. We discuss relationships between 
such Hamiltonians and the symmetry classes established for fermionic systems. We examine the 
density p{lo) of excitation frequencies uj, showing how the universal behavior piuj) ~ a;* for small 
Lo can be obtained both from general arguments and by detailed calculations for one-dimensional 
models. 
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An understanding of universal properties of excitations 
in disordered systems occupies a central place in con- 
densed matter physics. Much of the work in this area 
has focussed on systems of non-interacting fermions, as 
models for quasiparticles in disordered conductors or su- 
perconductors with interactions treated at the mean-field 
level. Classification of such systems according to symme- 
try provides an important starting point, and in many 
instances the possibilities are the three represented by 
the Wigner-Dyson random matrix ensembles A re- 
cent development, however, has been the recognition that 
there exist symmetry classes additional to those of the 
Wigner-Dyson ensembles. These additional symmetry 
classes arise in fermionic systems which have either chiral 
symmetry, as for tight-binding models on bipartite lat- 
tices with only off-diagonal disorder ||], or particle-hole 
symmetry, as for the Bogoliubov - de Gennes Hamilto- 
nian in disordered superconductors A characteristic 
feature of both cases is that cigcnstates at positive and 
negative energies are related in pairs and zero energy 
emerges as a special point in the spectrum. 

It is natural to anticipate that similar mathemati- 
cal structure may be important for the theory of non- 
interacting bosonic excitations or classical harmonic 
modes in disordered systems. Our aim in the following 
is to examine how far this is the case and what conse- 
quences it has. Studies of systems of this type have a 
long history, with celebrated early work by Dyson Q 
on the dynamics of a disordered chain of masses and 
springs, and applications which include phonons in dis- 
ordered solids and spin waves in random magnets . 
An obvious parallel between between excitations in these 
systems and those in fermionic systems belonging to the 
one of the additional symmetry classes is that bosonic ex- 
citations arise in positive and negative frequency pairs, 
and zero is a special point in the frequency spectrum. 
A second, more formal parallel is that random magnets, 
for example, in common with superconductors, may give 
rise to quadratic Hamiltonians containing terms that an- 
nihilate and create particle pairs. There are also clear 
differences between systems with bosonic and fermionic 
excitations, which prevent an elementary transcription 



of established ideas. Most importantly, while the exclu- 
sion principle guarantees that a fermionic system with 
a quadratic Hamiltonian has a ground state, for a cor- 
responding bosonic system constraints must be imposed 
to ensure that it is stable. As a consequence, even to 
construct a phenomenological treatment such as random 
matrix theory for bosonic excitations, it is necessary to 
keep in mind their origins in an appropriate non-linear 
problem. In turn, a distinction arises that is specific 
to bosonic excitations, between those which are Gold- 
stone modes and those which are not. We focus below 
on the latter and discuss elsewhere systems with contin- 
uous symmetry [Q. 

Any quadratic bosonic Hamiltonian can be written in 
the forms 

N 

Here, Xi and pi are the coordinates and momenta of the 
oscillators, = pi, S.N+i = Xi with 1 < i < iV, and 
the vectors at, a have as entries bosonic creation and 
annihilation operators aj , fli = (xiztipi) / \/2. The matrix 
Ti. is real symmetric, F is Hermitian and A is symmetric. 
The condition for time-reversal invariance is that C = 
or, equivalently, that F and A are real. 

We are concerned with spectral properties of models 
of this kind which are generic when Ti. is random. One 
anticipates that these will be found at small excitation 
frequencies uji, and we concentrate particularly on their 
average density p{lu) = N^^{Y^ 5(uj — uji)) at small uj. In- 
deed, the final expression of Eq. ([^) illustrates the parallel 
between bosonic models and the Bogoliubov - de Gennes 
Hamiltonian, from which universal behaviour has been 
derived for the density of quasiparticle states in disor- 
dered, gapless superconductors In fermionic systems 
such behaviour arises essentially from the interplay of 
disorder and symmetry, and does not depend on details 
of the distribution of Hamiltonian matrix elements. By 
contrast, for the bosonic systems we discuss, we show 
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that the requirements of stabihty impose features on the 
distribution of 7i which determine the form of p(w) at 
smah Lo. 

A simple argument leading to this conclusion has been 
given previously [p|-p^. In outline it is as follows. Recall 
that the Hamiltonian of Eq. (|l|) is characterised not only 
by the oscillator frequencies tOi but also by the eigenval- 
ues K„ of the matrix 7i, which we refer to as oscillator 
stiffnesses] let d{K) — {'^N)~^ (^5[k, ~ Kn)) be their av- 
erage density. From a discussion of the curvature distri- 
bution at absolute minima of random functions of one 
variable, it is suggested Q that ^(k) cx k^/^ for small 
K. Then, using the relation ujf — K1K2 which holds in a 
single- mode system, one arrives at p{uj) cx d{uj'^)u! cx 
for small ui. A weakness at both steps in this argument 
is that one degree of freedom is treated in isolation. We 
show here how the same behaviour emerges without such 
a restriction. 

In a stable system, stiffnesses are positive and frequen- 
cies are real. It is helpful to introduce a description 
that guarantees this property. To this end, write H as a 
square of a real matrix Q, in the form TC = Q^'^Q (pos- 
sible provided all En > 0). Also, note that frequencies 
are the eigenvalues of an auxilliary matrix Ti.' — a2'H, 
where (T2 = cry ® In, <^y is the usual Pauli matrix and 
Im is the N X N identity matrix |ll|]. Since the eigen- 
values of (J2Q^Q coincide with those of = Q(y2Q^ , 
which is Hermitian and antisymmetric, frequencies are 
real and occur in pairs ±tJi. While there is in general no 
simple relation between stiffnesses and frequencies, sev- 
eral special cases represent important exceptions. If (in 
Eq. (|^)) M = /jv and C = (as for vibrational problems 
with all masses equal), then = ±-y/K7, where the non- 
trivial stiffnesses appearing here are the eigenvalues of 
K\ whereas if M = K and C = (as for a random-bond 
ferromagnet) , then uoi — 

Moving beyond these special straightforward 
approach to multimode problems comes from concentrat- 
ing on disorder realisations in which one stiffness, say ki, 
is much smaller than all others. The probability den- 
sity for such disorder realisations varies as k^^^, by the 
arguments of Ref |8j. In the limiting case ki = 0, two 
eigenvalues of Vl vanish (since det(7i) = det(f7), and since 
the eigenvalues of f2 are paired), while for < ki <C k„, 
n ^ \, perturbation theory yields a pair of excitation fre- 
quencies uji cx ±y^. As a result, p(w) cx for a; <C Wp, 
where LOp is the lowest excitation frequency in a typical 
sample. More generally, consider a macroscopic system in 
which low frequency modes are localized with a localisa- 
tion length which remains finite as w — > 0. Treating each 
localization volume independently, we obtain p{lij) cx uj'^ 
for w <C Wp, where uop is in this case the lowest excitation 
frequency in a typical localization volume. 

Turning to detailed calculations, it is useful to con- 
struct a formulation in which Q appears linearly, by con- 



sidering the enlarged matrices 
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Clearly, the eigenvalues of Ti, and H' are squares of those 
of 7i and H' respectively. The 2x2 structure of ft dis- 
played in Eq(||) is interesting partly because it is the 
defining feature of the chiral symmetry classes, studied 
previously in a variety of contexts P,p^~p^ , while ma- 
trices with the structure of TL' constitute a new, chiral 
bosonic problem. Whereas much past work on chiral 
symmetry classes has taken the elements of Q to be inde- 
pendent random variables, drawn from a given distribu- 
tion, it is clear following our introductory discussion that 
a central concern in our case is to determine the distribu- 
tion of and hence that of Q. Moreover, to allow further 
progress, the form of Q must be sufficiently simple (for 
example, local), and it does not seem obvious in advance 
whether this will be the case for any given problem. In 
fact, there is considerable flexibility, since 7i determines 
Q only up to left multiplication by an orthogonal matrix, 
and we have found convenient, explicit expressions for Q 
in several one-dimensional examples In addition, the 
approach proves useful even for problems in which Q is 
known only implicitly, as we now show. 

We consider the random field XY model, which has 
been studied extensively as a description of pinned charge 
density waves ITtI . It has the Hamiltonian 



H = 



dx 
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where 11 is the momentum conjugate to the angle 0, 
h{4>,x) = h{x) cos{(j){x) — x(^))j s-nd h(x) and x{'^) 
the random amplitude and phase of an applied field, 
with {h{(j),x)) = and {h{(j), x)h{(t>' , x')) = S{x - 
x')hQ co8{(f) — (j)'). This continuum problem has, in the 
discrete notation used above, M = I^, C = and hence 
Lu = ±yfKi. We are therefore concerned only with K: 
for the sake of an obvious analogy, we refer to it as a 
Hamiltonian, denoting it by H and its eigenvalues by E 
in place of n. 

Stationary configurations of the field cj) satisfy 



+ (9^/1(0, a;) = 



(4) 



while amplitudes ip of normal mode excitations about the 
ground-state 0o, with frequency uj = ±\/E, obey 



Hijj = -dliP + dlh{(j)o, x)i: = EiP . 



(5) 



We want to find the density of states for the Hamiltonian 
of Eq. (^) . In this form, this problem has been the sub- 
ject of many publications p^^-pj[|. Summarising what is 
known: a characteristic energy scale, the pinning energy 
pO| Ep, separates two regimes. For E ^ Ep the den- 
sity of states is only weakly affected by the random field 
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and can be computed perturbatively in the field strength. 
For E <^ Ep the density of states is strongly influenced 
by disorder and believed to vary with i? as a power law, 
d{E) (X Ef^, but the value of (3 remains controversial 
|p7|,pp^. Here we show, in argeement with arguments 
due to Aleiner and Ruzin ^ and to Fogler that 
/3 = 3/2. 

Since Eq. (^) has the form of a Schrodinger equation 
with a random potential d'^h{<po,x), it is tempting to 
anticipate behaviour familiar from simple choices of po- 
tential distribution, such as Gaussian white noise. That 
would be too naive, however, because - as for any bosonic 
Hamiltonian - the spectrum of Ti. is positive and so the 
potential cannot be arbitrarily random. To investigate 
the implications of this fact, it is useful to follow Feigel- 
man ||l8|| , introducing the notion of a partial energy, 
£{(j),y)- the ground-state energy of the half-chain with 
coordinate values < x < y and boundary condition 
4>{y) = </*■ Interpeting Eq. as an equation of motion 
for a system with space coordinate and time coordi- 
nate X, £{(j),y) plays the role of an action and satisfies 
the Hamilton-Jacobi equation 

d,£ + ]^{d^£f ^h{^,x) (6) 

while the classical trajectory (t>Q{x) obeys dipo/dx — d^£. 
Now define V{x) = d'^£{(j)o{x), x), which we call the chi- 
ral potential. As a direct consequence of Eq. (^) it satis- 
fies 

A.V + V' = dlh{M^),x). (7) 

The Hamiltonian of Eq. ^ can therefore be rewritten as 
H = Q'^Q, where 

Q^-± + V{x). (8) 
dx 

Thus the bosonic problem of Eq. is equivalent to a 
one-dimensional chiral problem specified by Eqs. (^) and 
in which V{x) should be determined along with 
(j)o{x). 

At this point we can draw on the work of Comtet et 
al p^ , in which one-dimensional chiral problems of this 
kind have been analyzed in detail. In particular, they 
show that if the chiral potential V{x) has a positive av- 

1/2 

erage ( V{x) ) — Ep (sometimes referred to as stagger- 
ing Q), then the low-lying {E Ep) states of Eq. (||) 

— 1/2 

are localized with a localization length ^ ^ Ep . In 
this regime, the integrated density of states N{E) — 
d{E')dE' can be calculated as the probability of a 
negative fluctuation of V{x) in the interval xi < x < X2 
for which 2 dx V(x) = 2U < \og{E). The probabihty 
for a rare event of this kind is expected under many cir- 
cumstances to vary with E as exp[a log(i?)] = £'", but 
the value of a depends on the details of the distribution 



of V{x). In the following, we establish the equivalence 
implied by our notation, between ( V{x) ) and pinning 
energy, and calculate the low-lying density of states for 
the random field XY model by studying rare negative 
fluctuations of the particular chiral potential that arises 
in this context. 

To appreciate why ( V{x) ) > 0, consider the ground 
state of the half-chain as a function of the boundary con- 
dition, (p: it varies smoothly except for jumps at a small 
number of isolated points within one period. At these 
points £(0, y) is continuous but has an unward cusp [ p^ ; 
since £{4>,y) is periodic, an average of V{x) = com- 
puted excluding these cusps (as on (po) is naturally posi- 
tive, and its size can be estimated from the average sum 
of discontinuities in d^£. Alternatively, one can attempt 
a more detailed analysis of Eq. (|^), which is similar to 
the KPZ equation ||21|], but with two differences. One, 
the absence of a diffusion term Ddl£, is unimportant 
since it is well known that the limit D ^ in the KPZ 
equation yields solutions which are the global minimum 
of Eq. (H). The second, involving the nature of noise cor- 
relations is more significant: while studies of the KPZ 
equation deal with noise that is uncorrelated in both x 
and (j), we are concerned with noise that has long-range 
correlations in cj). This situation is familiar in the con- 
text of Burgers turbulence ||2^, and indeed, introducing 
u = dcf,£ one finds the Burgers equation Ux + uu^ — d^h. 
Cusps in 5(0, y) as a function of (j) corrspond to the shock 
waves of Burgers turbulence, and using dimensional anal- 

1 /3 

ysis the dependence on noise strength ( V{x) ) ^ can 
be obtained from Eq. (^, as in the standard picture of 
the pinning energy. 

Turning to fluctuations in V{x), a direct attack, via 
a solution of Eq. (^) , is not appropriate since many as- 
pects of Burgers turbulence remain controversial. Instead 
we treat Eq. (^ as a Langevin equation for V{x) with 
d'^h{(t)o{x), x) playing the role of a random force. If (j)o{x) 
is determined by minimising Eq. (|^) for fixed 0(0) and 
(j){L), the resulting random force has built-in correlations 
which ensure that V{x) remains bounded. Alternatively, 
one may solve Eq. (|^) for (t>{x) with fixed initial condi- 
tions, 0(0) and 9a;0(O). In this case, the random force is 
uncorrelated and appropriate values for dx4>{Q) generate 
ground states corresponding to different choices of (t>{L). 
Other values of 9a;0(O), however, correspond to maxima 
of Eq. (|^) , and still others to higher energy local minima. 
From the perspective of the Langevin equation, a max- 
imum of Eq. (H) is signalled by escape of V{x) towards 
large negative values. Those noise realisations that result 
in escape should be eliminated by supplementing Eq. (|^) 
with an absorbing boundary condition at V{x) = — oo. 
The trajectories of V{x) that are not absorbed should be 
weighted in two ways. First, we must eliminate solutions 
to Eq. (^ which are only local minima ofEq. (|). They 
are separated from the absolute minimum by maxima 
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and hence have neighboring trajectories on which V{x) 
escapes. The necessary weight is therefore the probability 
that no near neighbors of a given trajectory are absorbed. 
Second, we weight trajectories representing the absolute 
minimum for different dx(f>{0) in such a way as to give a 
uniform density on (l)(L). 

Omitting in the first instance both these weights, the 
form of the probability distribution P{U) at large nega- 
tive U can be found by passing from the Langevin equa- 
tion to a path integral, fixing U with a Lagrange multi- 
plier, and calculating the lowest eigenvalue of the corre- 
sponding Fokker-Planck operator, which is non-zero be- 
cause of the absorbing boundary condition. In this way 



we obtain P{U) 



^2U 



as U 



-oo. To calculate the 



weighting factors, we consider a family of trajectories pa- 
rameterized by s, with coordinate (j)(x)+r](s, x) and chiral 
potential V{x) + W{s,x), where ri{s,x) and W{s,x) are 
assumed small. Deriving linearized evolution equations 
by expanding Eqs. (Q) and (0), and solving these, we find 
■q{s,X2 



•q(s,xi) and W{s,X2) = e '^^F{s), where 



F{s) = W{s,xi) 

%xp(^3^ V{y)dy^dlh{(l),x)dx. (9) 

In order that no trajectories neighboring 4>{x) are ab- 
sorbed, we require W{s,X2) ^ —1 for all s and hence, 
at large negative U, F{s) > —e^^, which occurs with 
probability . Of the surviving trajectories, we should 
retain only those that have a fluctuation in the inte- 
grated chiral potential similar to U. These satisfy the 
condition W{s,X2) ^ 1, and hence F{s) < e'^^ , imply- 
ing |77(s,xi)| < and |77(s,X2)| < e'^^ , thus generating 
a weight e^^ . Combining these factors, the probability 
density that U is large and negative on 4>oix) varies as 
e^'^ , yielding N{E) ^ E^^'^, and hence p{uj) ~ w'*, as an- 
ticipated above. A detailed account of these arguments 
will be presented elesewhere 0. 

The mapping described here of Eq. (^) into a chiral 
problem also goes through for a discrete version of the 
XY spin chain. It gives, not surprisingly, a staggered 
one-dimensional random hopping model, as studied for 
example in Refs. This formulation is especially 

well suited to numerical calculations using the transfer 
matrix technique. Fig. 1 shows the result of a compu- 
tation of N{E) for a XY spin chain of 10^ sites, with 
Hamiltonian H = — cos(0i — where 

= Aicos(0 — Xi)i ^-iid [^cos(x), ^sin(x)] is uni- 
formly distributed on a disc of radius 0.01. The power 
law obtained analyticially is confirmed. 

We remark finally that, while the random field XY 
chain represents a special case in the sense noted above, 
since frequencies and stiffnesses are simply related, we 
expect the same form more generally for p{uj), as implied 
by the discussion preceeding Eq. ^ . Calculations for the 
random field Heisenberg spin chain lead to two-channel 



LU 




FIG. 1. The integrated density of stiffnesses, compared 
with the theoretically predicted power, 13+1 — 5/2 

chiral problems involving Ti. and Ti.' . For such problems 
we obtain 0| d{K) ^ k^^"^ and p{uj) ^ uj^. 
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